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ABSTRACT 


The  scattering  of  electromagnetic  waves  from  the 
ocean  surface  when  a  current  is  present  is  studied.  The 
electromagnetic  waves  are  singly  scattered  by  two-dimensional 
linear  surface  gravity  waves  interacting  with  a  surface 
current.  The  surface  is  described  in  terms  of  a  Wigner 
function,  i.e.,  the  spectral  density,  which  has  finite 
spatial  gradients.  The  effect  of  two  dimensions  on  the 
spectral  perturbation  is  explored  in  a  preliminary  calcu¬ 
lation  for  a  Phillips'  ambient  spectrum. 
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1 .  Introduction 


In  a  previous  paper  we  have  reviewed  the  effect  of 

surface  currents  in  modulating  the  ocean  surface  wave 

spectrum"*'.  In  this  work  we  have  restricted  ourselves  to 

the  linear  approximation  for  the  surface  waves.  The  role 

of  nonlinear  surface  wave-surface  wave  interactions  has 

2 

been  studied  in  a  preliminary  way  elsewhere  ,  and  a  code 
has  been  developed  to  make  a  more  thorough  study  of  non¬ 
linear  interaction. 

In  this  paper  we  use  the  rather  simple  results  of 
III  to  survey  qualitative  features  of  electromagnetic  wave 
scattering  from  the  surface.  A  simple  version  of 
scattering  theory  will  be  used  in  which  the  electro¬ 
magnetic  wave  is  scattered  only  once.  The  more  quanti¬ 
tative  aspects  of  the  scattering  will  be  addressed  in  a 
later  paper. 


2.  Description  of  the  Scattering 

Following  Callen  and  Dashen3  we  consider  the 
scattering  of  a  scalar  wave  field  by  ocean  surface  waves. 
They  separated  the  ocean  surface  into  a  long  and  short 
wavelength  part;  the  long  wavelength  part  having  a  small 
curvature  at  the  ocean  surface.  The  scattering  is  per¬ 
formed  by  the  short  wavelength  surface  gravity  waves 
riding  on  the  longer  waves.  In  the  present  analysis,  the 
long  wavelength  waves  become  the  surface  current.  This 
will  be  discussed  in  the  next  section.  For  scattering  an 
incident  plane  wave  of  wavenumber  k  to  a  final  wavenumber 
Jk* ,  we  have  the  scattered  field  in  a  simple  single  scattering 
approximation3 
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sc 
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5 (x»t) 


(2.1) 


Here  &k  =  k  -  k*  is  the  change  of  wavenumber  vector,  R0  is 

^  ^ 

✓ 

an  amplitude  factor,  and  the  integral  is  taken  over  the 
illuminated  area  AQ  of  the  ocean.  The  vector  x  lies  in  the 
horizontal  plane  of  the  ocean  surface  when  no  waves  are 
present  and  S(x,t)  is  the  vertical  displacment  of  the 
surface  at  x  and  time  t  due  to  wave  motion. 

We  use  the  notation  of  III  and  write 


(2.2) 


(2.3) 


Equation  (2.3)  provides  a  Fourier  representation  of  Z  in  an 

assumed  large  area  A  of  the  ocean  and  the  Fourier  coefficients 

o 

a(k)  are  functions  of  time. 

As  is  customary  in  changing  from  "box"  or  discrete 
to  "continuum"  normalization,  we  may  replace  the  sum  over 
discrete  k  by  an  integral  over  continuous  k.  with  the  sub¬ 
stitution 


V  Ao 


d2  k 


(2.4) 


The  power  spectrum  of  ?  (jx,t)  at  a  position  x  is 
4 

given  by  the  Wigner  transformation 


where  <•*•>  represents  an  ensemble  average  over  many 
observations  of  the  sea  state. 


(2.5) 
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<a(k)  a(p)>  =  <Z(x)  Z(y)>  =  0  (2.9) 

i . e . ,  the  phase  of  the  complex  amplitude  varies  so  rapidly 
in  time  that  its  net  contribution  to  the  statistical 
average  vanishes. 

If  we  wish  to  use  continuum  wavenumber  variables, 
we  refer  to  (2.4)  and  define  the  spectral  function 


(2.11) 


Now,  we  return  to  Eq.  (2.1)  and  obtain  for  the  mean 


scattered  power  the  expression 
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-  ( 2tt  )  2  jR°-L  /  d2x[y  (x,Ak)  +  'F  (x ,  -  %)]  •  (2.12) 


When  the  illuminated  area  AQ  is  sufficiently  small  that 
y  is  nearly  constant  over  that,  area,  we  have 

P  S  ( 2tt ) 2  Aq  I  R0 1  2 U5o  'A£>  +  ^o,  “  A&]  (2*13) 
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— ***< 


where  lies  within  the  illuminated  area.  To  observe 
¥  Ajc)  or  ¥£<o  -  Ajc)  individually,  it  is  necessary  to 

measure  the  Doppler  shift.  For  backscatter,  for  example, 
a  downshift  in  frequency  will  be  associated  with  ¥  (xq  Ajc) 

and  a  shift  up  for  ¥  (xQ  -  Ajc) . 

It  is  seen  from  Eqs.  (2.12)  and  (2.13)  that  the  • 
basic  requirement  for  evaluating  the  scattered  energy  from 
the  ocean  surface  requires  the  expressions  for  Z  obtained 
in  III  to  calculate  ¥. 


In  III  we  assumed  the  interaction  of  the  surface 
gravity  waves  with  the  prescribed  surface  current  to  be 
"turned  on"  at  time  t  =  0.  The  expression  obtained  for 
Z  was  of  the  form 


Eik-x 

e~~A(k)  Gk  (x,t)  ,  (3.1) 

Jc 

where  Gk(x,o)  =  1.  The  quantities  A(k)  thus  correspond  to 
the  a(k)  of  Eq.  (2.3)  when  there  is  no  surface  current. 

The  "ambient"  power  spectrum,  corresponding  to  the  absence 
of  a  surface  current,  is  then 


or 


A 

w  (X,k)  =  -2-*.  F  (x,k) 
a~~  (2u ) 2  a~~ 


(3.2) 


(3.3) 


When  the  ambient  wave  spectrum  has  a  negligible  spatial 
variation,  we  have 


Fa(^5S}  E  Fa(iS>  =  %<|A(k)r> 


(3.4 


Continuing  to  follow  III,  we  suppose  the  surface 
current  to  be  sinusoidal,  with  wavenumber  K  and  phase 
velocity  c ,  and  to  have  the  form 

U(£)  =  1  U(£)  =  1  UQ  cos(K£)  ,  (3.5) 

where 

5  =  x  -  Cj  t  .  (3.6) 

The  modulation  function  Gk  was  modelled  for  the 
"resonance  region"  in  III  as  follows.  First,  the  resonance 
wavenumber  kr  is  determined  by  the  condition 


ci  ~  cg(kr}  cos0r  =  0  '  (3.7) 

where  Cg(kr)  =  (g/k^y/2  and  cos0r=  Jcr»j./k.  The  resonance 

region  was  described  as  extending  over  a  range  jT 


Ak  =  K  S 


S  =  2 
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U 


CjO  cos  0-2) 
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(3.8) 


about  k._  . 
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A 

where  K  =  Ki  and  U  =  U  i.  The  relation 
—  — o  cr- 

00 

-  i 


iz  sin0 


(z)  , 


where  J^fz)  is  a  Bessel  function  of  the  first  kind  of 
order  v#  in  conjunction  with  Eq.  (2.5)  allows  us  to  write 
the  spectral  density  as 


(3.12) 
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Here  we  have  evaluated  G  at  the  wavenumber  Jc,  a  seemingly 
adequate  approximation.  The  quantity  Ji  in  Eq.  (4.3)  is 


f  ir- 

=  'Vis  [^(Cj  - 


Cg  COS0  } 


Uo  T  t  3  i  2r  1  ,, 

- 7  ct£  "  7  v  ca  cos  6  *  {4*4) 

;x  -  c  cos8  j  L  1  4  k>  g  J 


For  wavelengths  2r/k  much  less  than  the  wavelength 
at  which  resonance  can  occur,  there  is  another  contribution 
to  (4.3)  due  to  interaction  with  waves  of  wavenumber  I,  near 
resonance.  In  this  case,  Eq.  (4.3)  can  be  used,  but  with 

A 

G  given  by  Eq.  (5.15)  of  III: 


|g|2S  1  + 


Pcogfi.  +  fd2L  V  (L)  ^  (k»L 

2  Cj  J  a  ~  2  - — 


L)  Re(G-l)  ,  k  »  L.  (4.5) 


Here  we  have  assumed  that  the  ambient  spectral  function 
¥  (L)  is  independent  of  x,  as  was  done  in  III. 

a 


»  u 


5. 


Simplification  for  Small  Amplitude  Current  Effects 
At  this  point  it  is  convenient  to  make  the  trans¬ 
formation  (1.12)  to  the  continuum  power  spectra  yfe/k) 
in  Eqs.  (3. .\2)  and  (4.3).  We  shall  also  assume  that  the 
ambient  spectrum  is  independent  of  position,  so 


(x,k)  =  '?  (k) 


(5.1) 


When  -  |K|  «  |k|  and  the  current  is  sufficiently  weak,  we 
can  expand  both  the  ambient  spectra  and  Bessel  functions  in 
power  of  ~j.  On  keeping  only  the  lowest  order  and  linear 

term,  we  can  carry  out  the  sum  over  v  in  closed  form  with 
5 

the  relations 


00 

I'. 


(a)  =  1 


\)  =  -0C 


00 


(a)  =  a 


c’or  the  case  corresponding  to  Eg.  (3.12)  we  obtain 


loal  , 

(x,k)  =  Qc)  <  1  +  1  — 2  -  1  cos^0sin  (Kg) 

+  igCOo**-)2  cos(KC)  jl  -  —j2/S 
jt|4  sin[2S1/8  sin(K£)]  'l 


+  2  — 
1TP; 


ti4  sin[2S^8  sin  I 


smK£ 


-  [k  K  sin<K«>]  • 


(5.2) 


For  waves  far  from  resonance,  we  make  the  corres¬ 
ponding  approximation  in  Eq.  (4.3)  to  obtain 


¥(x,k) 


/ 

=  V#  1 

k. 


2 (CI~Cg  COS0) 


2  cI(cos0-2) 


+  2c  COS0(1-COS0)  +  c  cos  0 
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r  3Vin 

U  cos0  i 

L  “i^TJ 

c-c  COS0 

L  I  g  J 

(5.3) 


Here  c  -  (g/k)  */2  and  cos0  =  i*k,  as  before, 
g 

Finally,  for  wavelengths  very  short  compared  to 
wavelengths  near  resonance,  we  use  Eq.  (4.5)  to  obtain 

V(x,k)  =  V  (k)  1  +  —  (cos0-2) 

^  ~  a  —  ^  2CJ 

*  JAts!l}  ^  (£•£) 


(5.4) 
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Thus 

3^  (k) 

k  --f-r--  =  -4  cos0  ¥  (k)  .  (6.2) 

3kx  a  ~ 

Our  parameters  are  chosen  to  correspond  to  typical 
oceanographic  conditions: 

_2 

U  =  10  m/sec 

o 

-2  -1 
K  =  10  m 

k  =  8.66m 
r 

Cj  =  0.50  m/sec 

0  =  20°  (6.2) 

Here  kr  is  the  resonant  wavenumber  defined  by  c ^  =  c^  cos0. 
The  strength  parameter  S  is 


V 


S  =  4.609x10 


(6.3) 


t r  =  Z - o~  =  23.09  sec  . 

res  o 


(6.4) 


Rosenbluth's  "pile  up  time"  is  then 


T  =  3894  sec 
P 


(6.5) 


and  the  extent  of  the  resonance  region  ^Eq.  (3.8)J  is 


Ak  =  KS*5  =  2.147  m' 


(6.6) 


For  the  spectrum  (6.1)  the  term  involving  an  integral 
over  £  in  Eq.  (5.4)  is  negligible  compared  with  the  other 

two.  In  Eqs.  (5.2),  (5.3)  and  (5.4)  the  term  involving 

r  7 

k  appears  to  dominate  ^except  for  t  >>  Tr  in  Eq.  (5.2)1. 

X  7 

This  was  pointed  out  by  Milder  in  his  discussion  of  spectral 
perturbations.  For  a  spectrum  of  the  form  k^  the  relative 
contribution  of  the  first  and  last  terms  in  Eq.  (5.4)  is 
approximately  —  ,  which  in  two  dimensions  is  .125. 

For  the  parameters  given  by  Eqs.  (6.2)  -  (6.6), 

Eq.  (5.7)  becomes,  for  an  adverse  current, 


f  2  2 

[k.)  |  1  -  .01  sinK£  -  4.34xlo~5(~)  cosK?  ■ 

V.  r  r  P  , 


(6.7) 
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The  term  linear  in  t  above  results  primarily  from  the 

spectral  gradient  term,  i.e.,  90%  of  its  value. 

To  illustrate  the  spectral  perturbation  in  the 

regions  away  from  resonance,  we  graph  Eg.  (5.3)  in 

Figure  1  excluding  the  interval  Ak  about  k^  .  It  is 

evident  that  abuve  the  maxima  of  the  surface  current 

(C  =  0) ,  the  longer  waves  are  enhanced  and  the  shorter 

waves  suppressed.  The  magnitude  of  the  perturbation  is 

found  to  increase  as  one  goes  to  smaller  angles.  At 

0  =0°,  the  perturbation  extremes  are  in  substantial 

8  9 

agreement  with  previous  calculations  '  . 

The  width  of  the  resonant  region  in  Figure  1  is 
given  by 

Ak  =  KS3* 

so  that  using  Eq.  (3.8)  with  the  parameters  we  have  been 
discussing, 

JT-  =  ■  ,-11-5— -  (6.9) 

r  >/cos20-.666 

which  is  plotted  in  Figure  2.  We  see  from  the  slow  rise 
of  Eq.  (6.9)  from  its  minimum  at  0  =  0°  that  the  resonance 
width  is  proportional  to  the  resonance  wavenumber  for 
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Figure  1, 
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-1.2 


|0|  -  30°  to  within  a  factor  of  2.  In  fact,  it  is  only 
within  a  relatively  narrow  angular  range  about 
0  =  +  35.26°  that  the  two-dimensional  effects  become  very 

strong . 

Finally  we  assume  k  >>  kres  and  write  Eq.  (5.4)  as 

'l!(x,Jc)  3  S'  (Jc)  [l  -  0.06  cos  (KS)  ]  •  (6.10) 

Here  also,  the  3<T /3k  term  provides  the  dominant  contri- 
bution  to  the  modulation.  The  coefficient  (.065)  in 
Eq.  (6.10)  is  seen  to  be  the  asymptote  of  the  relative  per¬ 
turbation  for  the  short  waves  as  seen  in  Figure  1. 

in  Figure  1,  we  see  the  effect  of  the  dependence  of 
the  resonant  wavenumber  (kr)  on  angle,  and  from  Eq.  (6.9) 
the  dependence  of  the  spectral  interval  on  angle.  The 
solid  line  is  the  spectral  perturbation  for  9  =  20°  and 
the  dashed  for  6=0°.  The  resonance  interval  (kr  ±  Ak/2) 
is  indicated  for  the  8  =  20°  case.  The  0  perturbation 
is  seen  to  be  shifted  to  a  higher  wavenumber  and  the 
maximum  perturbation  exceeds  the  20°  case,  although  for 
k  —  (kr  -  Ak/2) g_20°  the  20°  Perturbation  exceeds  the  0° 
perturbation.  The  suppression  of  the  shorter  waves,  how- 
ever,  is  greater  for  all  wavenumbers  m  the  0  case. 
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